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We provide a scheme for the generation of entangled number states of Bose-Einstein condensates 
in multiple wells with cyclic pairwise connectivity The condensate ground state in a multiple 
well trap can self-evolve, when phase engineered with specific initial phase differences between the 
neighboring wells, to a macroscopic superposition state with controllable entanglement - to multiple 
well generalization of double well NOON states. We demonstrate through numerical simulations 
the creation of entangled states in three and four wells and then explore the creation of "larger" 
entangled states where there are either a larger number of particles in each well or a larger number of 
wells. The type of entanglement produced as the particle numbers, or interaction strength, increases 
changes in a novel and initially unexpected manner. 
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I. INTRODUCTION 

Entanglement, a nonclassical correlation between two 
or more physical systems, lies at the heart of the profound 
difference between quantum mechanics and a local clas- 
sical description of the world [1]. Apart from their dis- 
cussions in the philosophical and foundational aspects of 
quantum mechanics [2], entangled states in recent years 
have become an essential resource for the emerging field 
of quantum information processing [3]. Entangled and 
squeezed states hold promise in studies related to quan- 
tum measurement, Heisenberg limited atom interferom- 
etry and precision measurements [4, 5], and quantum 
computing and quantum communication. Since multi- 
particle entangled states can be more useful than the 
two-particle/two-photon states, there have been steady 
attempts toward creating such states [6]. Such states 
have been created with several systems - with five pho- 
tons [7] , eight atoms in an ion trap [8] , ten nuclear spins 
in a molecule [9], and also with cold atoms in an optical 
lattice [10, 11]. There have been several proposals to cre- 
ate superposition states with a large number of particles 
such as using a tiny mirror [12] and microorganisms such 
as viruses [13]. 

While the consequence of entanglement for an 
Einstein-Podolsky-Rosen (EPR) pair is quantified in 
Bell's inequality [14], a more striking conflict between 
quantum mechanics and local realism is exhibited by 
three maximally entangled spins also known as the 
Greenberger-Horne-Zeilinger (GHZ) states [15]. GHZ 
state of N spins has the form 
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This state can also be written in the notation |00...0) + 
1 11... 1) where |0) and |1) are the basis states. The su- 
perposition of two macroscopically distinct states, rather 
than simply the internal degrees of freedom, each occu- 
pied by all TV particles, has been discussed by Schrodinger 



in the famous cat parable [16]; partial realization of such 
Schrodinger's cat states has been obtained with Joseph- 
son junction loops [17, 18]. That macroscopic superpo- 
sition states are highly entangled has been discussed in 
Ref. [19]. 

The analog for the GHZ state for TV particles in two 
wells is denoted like Eq. 1, -^(|JV,0) + |0,TV)), and is 
colloquially referred to as "NOON" states. A general- 
ization of this two state model to multi dimensions has 
been discussed in Ref. [20]. In this paper, we discuss 
the generation of macroscopic entangled number states 
of a multiwell Bose-Einstein condensates (BEC) of the 
approximate form 
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This is the multiwell generalization of the double well 
NOON states where a macroscopic number of particles 
are simultaneously in M different locations, with M > 2. 
BEC in optical lattices [21] has been a promising research 
area with many new observations such as the superfluid 
to Mott insulator transition [22] and number-squeezed 
states [23]. Superfluid and Mott insulator states are the 
ground states of bosons in an optical lattice, whereas the 
multi-positional Schrodinger cat state of type Eq. 2 is 
the highest lying excited state. Due to the coherence 
properties and versatility of cold atom systems, it may 
be an ideal system to create such entangled states. 

We show that states approximating the extreme entan- 
gled states of Eq. 2 may be generated in a controlled fash- 
ion by time evolution of appropriately phase imprinted 
ground states of a multiwell BEC with periodic bound- 
ary conditions for M = 3,4, and 8. The physical mecha- 
nism for creating such states can be understood from the 
phase space picture of a double well BEC [24, 25] where 
it was shown that the ground state wave packet displaced 
in phase and put on a hyperbolic fixed point of its un- 
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derlying semiclassical phase space dynamically bifurcates 
to a macroscopic superposition state, a highly entangled 
state. Similarly, for a multi-well BEC, phase imprinting 
the ground state moves it to an unstable equilibrium, and 
subsequent dynamics creates entangled states. We show 
that the choice of initial barrier heights, which deter- 
mine the extent of ground state number squeezing, and 
the rate of barrier ramping can be used to control the 
entanglement of the final states. 

Based on results obtained for two, three, and four 
well configurations, we conjecture a generalized formula, 
for M wells, for the phase offset between neighboring 
wells appropriate for the generation of number entangled 
states. Finally, we extend our analysis to larger systems 
where there are a larger number of particles in each well 
or a larger number of wells. In these cases, we find sur- 
prising results indicating the formation of a new type 
of number entangled state in four wells. We study the 
impact of increasing the number of particles or the inter- 
action strength on the resultant entangled state. 

The paper is organized as follows. In Sec II, we intro- 
duce our model and discuss the numerical methods. In 
Sec III, we review previously published [24] results for 
the generation of NOON like entangled states in a dou- 
ble well. This is done in order to motivate and generalize 
the double well results to multiple wells. In Sec IV, we 
present our study of three and four wells showing the 
phase engineering method for creating entangled states. 
We focus on small number of particles, mainly showing 
that the method works for multiple wells, and that we 
can control the final state by controlling the initial bar- 
rier height and barrier ramps. In Sec V, we study this 
for larger number of particles in four and eight wells, and 
find that a new type of entangled number state is also 
created in the process. Finally, we summarize our results 
in conclusion in Sec VI. 



II. METHODS 
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FIG. 1: Shown are multi-well configurations for three, four, 
and eight wells with periodic boundary conditions in a one 
dimensional circular array. 

The physical configuration we study assumes multiple 
wells connected in a circular array, which has been real- 
ized experimentally by [26], and is illustrated in Fig. 1. 

We approximate the physics of a BEC in a multi- 
well potential by the Bose-Hubbard (BH) model [27, 28]. 



Thus 
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where n* = a J a* is the number operator, J is the nearest 
neighbor tunneling term, U is the on-site energy, and e* 
is the energy offset of the ith lattice. To simplify a theo- 
retical study, we make a one parameter approximation of 
the tunneling and interaction strength: U/J = l/e~ a ^; 
and for the symmetric wells explored here, e* = 0. Be- 
cause of the tunability in optical lattices, U and J can 
be changed in time with varying lattice depth, and we 
take a(t) as a function of time t. a(t) is a dimensionless 
parameter that can be mapped onto the barrier height. 
This parametrization allows a simple study of continu- 
ous change of barrier height through the variation of a 
single parameter a(t). For example, for a lattice made of 
red detuned laser with A = 985 nm and for 23 Na, a bar- 
rier height 15E R gives U = 0.15E R and J = 0.07 E R [28] 
where Er = is the recoil energy from absorption of a 
photon; these experimental parameters then correspond 
to a = 2.14. 

Our numerical studies focus on three size regimes: I) a 
small number of particles in a small number of wells, II) a 
large number of particles in small number of wells, similar 
to the experiments of [23] , and III) a small number of par- 
ticles in a larger number of wells, such as the experiments 
of [22]. For TV identical bosons in M wells the number 
of Hilbert space dimension is D = ^^^rx^ • The latter 
two cases result in large D and necessitate the use of par- 
allel processing techniques. The largest systems we have 
investigated at the time of writing is 512 particles in four 
wells with 202 million nonzero entries, 12=22,632,705 and 
24 particles in eight wells, 12=2,629,575 with 3.5 million 
nonzero entries. Details of our parallel implementation 
of the Bose Hubbard model can be found in [29]. 

The main equation we solve is the time dependent 
Schrodinger equation (TDSE): 

Y}H) jk c k {t)=ih d ^- (4) 

k 

where H is the discretized BH Hamiltonian in the Fock 
state basis \m). The solution of the time dependent 
Schrodinger equation in the full Fock space is then of 
the form 

D-l 

\*(t))=Y, c m(t)\m) (5) 

m=0 

III. ENTANGLED NUMBER STATES OF THE 
BEC IN TWO WELLS 

In order to generalize the creation of entangled states 
in multiple wells, we briefly review here the physical prin- 
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FIG. 2: Comparison of the classical nonrigid physical pendu- 
lum phase space with quantum phase space shown with the 
Husimi probability distributions for double well eigenstates. 
Shown are (a) classical energy contour. Husimi projections for 
(b) ground state which is minimum uncertainty wavepacket 
centered at the origin, (c) 6th state, (d) 12th state, and (e) 
35th state which is analogous to a superposition of classi- 
cal pendulum rotor motions in two opposite directions. The 
analogy between the double well BEC and physical pendulum 
points to a way to the creation of macroscopic superposition 
states by displacing the ground state to the unstable equilib- 
rium points in the classical phase space. 

ciples behind creating such states in a double well which 
was described in Ref. [24, 25]. The physics of creating 
such states in a double well as well as its extension to 
multiple wells as presented here can be understood in 
terms of the underlying classical phase space dynamics. 

The most general state vector in a double well is a 
superposition of all the number states 

N 

l*>= E c{ n\W,N-n L ) (6) 

n L =0 



TV the total number of particles. Finding the eigenval- 
ues and eigenvectors of the Bose-Hubbard hamiltonian 
for two wells can be easily accomplished by diagonaliz- 
ing a (TV + 1) x (TV + 1) tridiagonal matrix, getting the 

coefficients Cn\. 

After we implement the cat state generation method, 
states of the following form can be generated, 

\*) = -L(\N-n L ,n L ) + \n L ,N-n L )) (7) 

When til = or TV, it becomes an extreme superposition 
state 

{^extreme) = ^= (\N, 0) + |0, TV)) (8) 

where all the particles are simultaneously in the left 
and right wells. Here the cat states are positional 
NOON states where TV particles occupy spatially sep- 
arated modes in two wells. There can also be NOON 
states of other kinds such as with TV particles occupying 
the two quasi-momentum modes of counter-propagating 
superfluid flow states in a rotating ring lattice that have 
been proposed in Ref. [30, 31]. 

Anderson [32] showed that the Hamiltonian for a sys- 
tem of two quantum fluids connected by a tunneling junc- 
tion can be described as a physical pendulum. For BEC 
in a lattice, in the semiclassical limit valid for large TV, 
the operators &i in Eq. 3 can be approximated by the 
c- numbers ^Jnle %di , where rii and 0i are the number and 
phase of particles in the ith well. For a two site Bose- 
Hubbard model, the Hamiltonian then turns into a clas- 
sical Hamiltonian of a nonrigid physical pendulum with 
the number and phase differences (n = (til — n^)/2,^ = 
@l — Or) between the wells as conjugate variables. The 
dynamics of double well BEC is then described by a clas- 
sical pendulum model that has been studied in detail in 
Ref. [33] and experimentally verified in [34]. Fig. 2(a) 
shows the classical energy contour showing the classical 
phase space structure. This system has two fixed points 
- (0,0) and (0,7r). The (0,0) is a stable equilibrium, while 
the (0,7r) is stable in the tt- state regime (UN/ J < 1) and 
unstable otherwise (UN/ J > 1). Two types of pendulum 
motions - oscillations and rotor motions appear in the 
phase space, below and above the separatrix. Question 
can be raised on how much of the classical phase space 
is actually contained in a full quantum analysis. This 
was answered in Ref. [24] showing quantum-classical cor- 
respondence between the classical pendulum and double 
well BEC. 

Husimi probability distribution [24, 35] can be used to 
project, in a squeezed coherent state representation, the 
classical phase space properties from quantum wavefunc- 
tions. In (n, 0) representation, Husimi function is defined 
as 



where is the number of particles in the left well, and 
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Here n' = n1, 2 Ufl , rather than being the simpler left par- 
ticle counter, and c n > is the corresponding Fock- state co- 
efficient. The 'coarse-graining' parameter k determines 
the relative resolution in phase space in the conjugate 
variables number and phase. 




FIG. 3: A visual depiction of the underlying physical princi- 
ples for the generation of entangled states in a double well. 
Quantum phase space pictures of Husimi projections are plot- 
ted here - (a) the ground state at t=0, (b) 7r-phase imprinted 
ground state at the hyperbolic fixed point, (c) at t=0.012 
ms the wavepacket is bifurcating along the separatrix, (d) at 
t=0.019 ms it continues to move along the separatrix, (e) with 
simultaneous ramping of the barrier the wavepacket splits 
completely at t=0.49 ms, and (f) at t=2.84 ms a sharply 
peaked entangled state is obtained which is a macroscopic 
superposition of particles simultaneously in both wells. 

Figs. 2(b)-(e) show representative Husimi projections 
for 40 particles for the ground state, 6th , 12th and 35th 
states respectively. The ground state is a wave packet 
centered on (0,0) in the classical phase space, with a fi- 
nite width in number and phase differences. The higher 
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FIG. 4: Shown is the evolution to a cat state in Fock space 
at the same time instants as in the previous figure: (a) the 
phase imprinted ground state at t=0, (b) at t=0.012 ms the 
Gaussian distribution broadens, (c) at t=0.019 ms it bifur- 
cates, (d) at t=0.49 ms it splits completely, (e) at t=2.84 ms 
a highly entangled state is formed; (f) is a magnified version 
of (e) showing the nonvanishing Fock state coefficients. 



lying states in (c) is a harmonic oscillator like state, (d) 
is a state which lies on the separatrix and (e) is a cat 
state, which in classical sense, is a superposition of clock- 
wise and counter-clockwise pendulum rotor states. After 
understanding this quantum-classical correspondence, we 
can argue that a ground state wave packet displaced in 
phase by tt and put onto the unstable fixed point (0,7r) 
would bifurcate along the separatrix and create a super- 
position of two pendulum rotor states, if allowed to time 
evolve. The tt displacement of the ground state can be 
accomplished by phase imprinting one of the wells by 
an amount tt which in experiments could be done by a 
phase engineering technique that has been demonstrated 
in experiments investigating solitons in the BEC [36]. 

For a concrete example of our method, we show results 
of a numerical simulation in Figs. 3 and 4. Fig. 3 demon- 
strates the underlying physical principles quite visually 
in quantum phase space using Husimi projections. First 
in panel (a) we have a ground state that is at the cen- 
ter of phase space, then in (b) we displace this by tt in 
the horizontal direction so that it is on an unstable equi- 
librium. Since it is no longer an eigenstate, it will time 
evolve spontaneously, and in this case the trajectories 
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in phase space follow along the separatrix and symmet- 
rically spits in two directions as in (c) at t=0.012 ms. 
In this process the phase space points reach the top and 
bottom of the separatrix as in (d) at t=0.019 ms. If we si- 
multaneously increase the barrier height during this time 
evolution, we can completely split the top and bottom as 
shown in (e) at t=0.49 ms, finally giving rise to a desired 
cat state in (f) at t=2.84 ms. The times for this exam- 
ple of double well as well as three wells in next section 
are given for a 87 Rb condensate, A = 840 nm, a sc = 5.8 
nm, J = 0.04£ , j Re _a where Er = is the recoil energy 
from absorption of a photon, and a is a dimensionless pa- 
rameter that can be mapped onto the inverse tunneling 
rate, and taking U = O.OAEr as approximately constant 
for calculation purposes. Here a varies as a = 3 + 2t. 
The evolution to an entangled state is shown in Fock 
space coefficients in Fig. 4 - the process of the formation 
of a superposition state can also be understood from this. 
For a two-component spinor Bose gas, Ref. [37] shows the 
generation of entangled state, where there is no need for 
initial phase imprint. We demonstrate here the genera- 
tion of number entangled cat states of the NOON type; 
we do not discuss generation of phase cat states [38] that 
are superpositions of many coherent phase states. Much 
of the intuition for later sections is derived from an in 
depth study of the double well as briefly described here 
and presented in detail in [24]. 



IV. ENTANGLED NUMBER STATES IN 
THREE AND FOUR WELLS WITH SMALL 
PARTICLE NUMBER 

In order to gain insight into the multiwell Bose- 
Hubbard model, we first analyze the quantum mechani- 
cal properties of the simplest multiwell potential, M = 3, 
assuming three symmetric wells in a circular array [39]. 

The state vector for three wells is a superposition of 
all the number states 

N 

l**>= E c&naWIm, na.ns) (U) 

n\ ,n2=0 

Here ni, r&2, and ns = N — ri\ — rii are the number of 
particles in each of the three wells. Fig. 5 shows the 

Fock space probabilities, Cn\,n 2 , for representative sta- 
tionary states for N = 12 and a = (U/J = 1). n\ 
and 77,2 are the Fock state indices and the vertical axis 
shows probabilities. The ground state in Fig. 5(a) is a 
broad Gaussian while the higher lying states, Figs. 5(c)- 
(d), are number entangled states of increasing extremity 
corresponding to increasing numbers of particles simul- 
taneously in all three wells, the highest of which in panel 
(d) approximates an extreme superposition state of the 
form |A/", 0,0) + |0,7V,0) + |0,0,7V). Note that there are 
still some nonvanishing Fock space components. As U / J 
increases the highest lying state approaches the extreme 



1 C nln2 1 




1 C nln2 1 2 




FIG. 5: (color online). Fock state coefficients for 12 particles 
in three wells: (a) the ground state, (b) 10th, (c) 76th and 
(d) 91st, the highest state. The ground state has a Gaus- 
sian shape, while higher lying states are entangled number 
states, m and 712 are the Fock state indices and the vertical 
axis shows probabilities. Points beyond the cross- diagonal are 
unphysical. 



state of Eq 2. The number of non vanishing Fock state 
coefficients determines sharpness, and thus (d) is sharper 
than (c). 

It is unlikely that such maximally entangled states can 
be generated via a sequence of single particle excitations. 
They may however, be dynamically generated via phase 
engineering from the appropriate ground state, as eluci- 
dated in previous section for the double well. Writing 
phases on part of a condensate is experimentally feasi- 
ble via interaction with a far off-resonance laser [36] , and 
is assumed to be sudden with respect to the dynamics 
of the condensate. Mathematically, this corresponds to 
multiplying the coefficients in an expansion of the type of 
Eq. 5 by e mi6>i , where |ni, r&2, ...rii, ...) is the correspond- 
ing Fock state, and 0i is the phase for particles in the ith 
well. 

Entangled state generation, obtained via integration of 
the linear time-dependent Schrodinger equation, is shown 
in Fig. 6, following phase imprinting of an initial phase 
difference of between the neighboring wells, and a si- 
multaneous linear ramping of the barrier as a = 0.5 + 1, 
as shown in Fig. 6(a) (t here is dimensionless). Panels 
6(b) shows the initial ground state; 6(c) at time 0.43 ms, 
the distribution broadens; 6(d) at 0.74 ms, in the pro- 
cess of splitting the state towards the three corners; and 
6(e) at 1.36 ms a sharp, although not extreme, entan- 
gled number state with its signature of three major non 
vanishing expansion coefficients. 

When an appropriately entangled state is reached the 
barrier is suddenly raised to halt further evolution in n- 
space. For the parameter values used here, a simple time 
evolution without any change of barrier also produces an 
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FIG. 6: (color online). Evolution to an entangled Fock space 
state for three wells: (a) barrier ramp showing the location of 
the following time evolved states: (b) initial state, (c) at 0.43 
ms the Gaussian distribution broadens, (d) at 0.74 ms the 
distribution is 'splitting', (e) A three-peaked state is formed 
at 1.36 ms; a macroscopic superposition of definite number 
of particles simultaneously in all three wells. A comparison 
with the double well Fock state evolution in Fig. 4 shows that 
the physical mechanism for generating entangled states is the 
same in a double well and multiple wells. 



entangled state, however barrier ramping is used here to 
sharpen the resulting state, and completely split Fock 
state coefficients into three parts. Control of the extrem- 
ity of the states can be achieved by choice of the ini- 
tial barrier height controlling the initial squeezing of the 
ground state. This is demonstrated in Fig. 7, where dif- 
ferent initial squeezing have been used for rows (1), (2) 
and (3). The columns show: (a) the ground state, and 
(b) the final state at the end of the barrier ramping. It 
is important to be able to tune to less extreme entangled 
states, as such states are more robust to loss and deco- 
herence [24, 25]. In our study here, we show the proof of 
principles that cat states can be produced and controlled. 
To be able to get extreme superpositions, analysis can be 
made with optimal control theory on the correct param- 
eters and ramping to be used. Phase imprinting with a 
phase difference of ^ produces an equivalent state, with 
different phase space dynamics. 

What is so special about phase imprinting ^? Sim- 
ilar to the simple double well, the triple well, M = 



FIG. 7: (color online). Entangled states evolved from ground 
states with different initial squeezing. Row (1) shows the 
states with a — 1.0 + 1: (a) initial ground state and (b) final 
state. Row (2) is for a — 0.5 + £ and (3) is for a — t. Column 
(b) gives the states at t=1.85 ms, t=1.36 ms, and t=0.99 ms 
respectively. The initial squeezing of the ground state thus 
determines the extremity of the resulting entangled states. 



3, can be thought of as two coupled pendulums [40] 
with complicated dynamics - quantum and semiclassi- 
cal aspects of three well BEC have been elucidated in 
Ref. [41-43]. Here the semiclassical conjugate variables 
are An i2 = n 1 - n 2 , An 2 3 =n 2 - n 3 , A6 12 = 1 - 2 and 
A#23 = #2 — #3- The unstable fixed points in these con- 
jugate variables are (0,0,^^,-^p 1 ), k = 1,2. So, a phase 
imprint of ^ or ^ puts the ground state wavepacket on 
the unstable equilibrium, and subsequent dynamics gives 
rise to these states. 

One important aspect of our method is the control- 
lability of the final state which works for multiple wells 
as illustrated in Fig. 7. It works in three ways - I) A 
simultaneous ramping of the barrier with the natural dy- 
namics at the unstable fixed point has been empirically 
found to be useful in directing the desired evolution of the 
wavepacket. II) Initial barrier height, that is the initial 
squeezing, helps shape the initial wave packet stretching 
it into different regions of accessible phase space; and, III) 
the initial barrier height sets the (negative) curvature of 
the potential at the hyperbolic fixed point, controlling 
the rate of splitting of the wave packet. 



7 





FIG. 8: (color online). Four well stationary and time evolved 
states: (a) ground state, (b) an entangled state evolved from 
the ground state following a tt relative phase shift. The three 
dimensions show the Fock state indices m, 712 and 713, prob- 
abilities are shown in the color intensity scale. For graphical 
clarity, only the points higher than 40% of the highest proba- 
bility are shown, with the highest probabilities normalized to 
1. 
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Gaussian. After a phase imprint of tt between neighbor- 
ing wells, the ground state can be evolved into an extreme 
entangled number state as shown in Fig. 8(b) for N = 16, 
U = 0.01 E R , J = 0.04 E R e~ a , with a = 0.175, and as- 
suming the case of 87 Rb in the previous example. Phase 
difference of tt between neighboring wells is equivalent to 
writing alternating tt phases on the lattice. The other 
fixed point dynamics of | and ^ also lead to symmetric 
states but do not lead to the highly entangled states of 
the kind we discuss in this article. 

In Ref. [44], a truncated Wigner approximation was 
used to study an alternating tt phase difference dynam- 
ics for even number of wells. In comparing our four well 
results to theirs, we have done an exact time evolution 
study, and find that the tt configuration in an even num- 
ber of wells that they have identified is just a special 
case of many phase imprint dynamics that could gener- 
ate interesting correlated states in multiple wells. Their 
changes in system parameters is to drive the system from 
stability to a regime of instability. On the other hand, 
we take our system to be in the unstable regime and 
demonstrate the controllability of entangled states with 
barrier manipulation; potentially useful for experimental 
detection. 

For the two, three, and four wells, we find M — 1 dis- 
tinct phase differences between the neighboring wells for 
the multi-well fixed points [40] . These are given by a gen- 
eral formula where j = 1,2,....M — 1, with M being 
the number of wells, which gives a tt phase difference for 
the M = 2 double well, a ^ and ^ phase difference 
for the M = 3 triple well, and a ^, 7r, and ^ phase 
difference for the M = 4 quadruple well configuration - 
we investigated the dynamics generated by all of these 
phase difference imprints. Note that the total change in 
phase in the circular loop is a multiple of 27r, a vortex 
like condition. We thus propose a general formula for M 
wells, 



(M-l) phase configurations for M wells 

FIG. 9: Initial phase configurations that generate number en- 
tangled states. These phase differences are motivated by the 
fixed point of the underlying semiclassical phase space - bifur- 
cation characteristics of the unstable equilibrium generating 
the macroscopic superpositions. Shown here are examples for 
a double well (a), triple well (b and c), and four wells (d, e 
and f). There are M — 1 different possibilities for M wells. 
The known case of tt phase differences for an even number of 
wells is seen to be a special case. 



We found that all the features of the double well en- 
tanglement generation apply to the three well case, and 
thus, many of the insights from the two and three well dy- 
namics can be extended to arbitrary number of wells in a 
circular array. Next, we explore it for four wells. For four 
wells, the values of phase imprints that take the ground 
state to fixed points are 77, tt and Fig. 8(a) shows a 
typical ground state in four wells that is an approximate 
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for the constant phase offset between neighboring wells 
leading to the dynamical generation of entangled states. 
Here j = 1, 2, .., M — 1, and Eq. (12), being valid for any 
number of wells, even or odd, provides a substantial gen- 
eralization of the tt phase offset mentioned in Ref. [44], 
which is valid only for the special cases of an even num- 
ber of wells and for j = M/2. Fig. 9 shows the phase 
configurations of Eq. (12). The multiplicity of Eq. (12) 
is prominent for large number of wells, e.g. for 12 wells, 
there are 11 phase offset possibilities. Symmetries may 
prevent all the imprinting offsets of Eq. 12 from generat- 
ing independent dynamics. 

Although detailed work on four and eight wells was 
done more recently and is being presented here as a co- 
herent whole for multiple wells, the work on three well 
was done much earlier as documented in Ref. [25, 45]. 
Many studies of quantum, semiclassical and entangle- 
ment aspects of three coupled BECs have appeared since 
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then, that continue to explore the rich yet simpler lattice 
physics of the system of three well BEC. 



V. LARGE ENTANGLED NUMBER STATES IN 
MULTIPLE WELLS 

In this section we investigate the creation of large en- 
tangled number states where there are a large number 
of particles in each well and/or a large number of wells. 
Following the techniques described above and in [29], we 
explored the creation of entangled number states in four 
and eight wells with periodic boundary conditions. We 
present results from our simulations of entangled number 
states in four wells in the Section V A, eight wells in V 
B, and in V C we analyze the types of entangled states 
generated and describe a new type of entangled number 
state. 

In systems with wells M > 4, visualization of multi- 
dimensional entangled states becomes more difficult. To 
facilitate the visualization we introduce the joint proba- 
bility function 

1 N 

P(np,n 7 ,t) = — K 1 ,n 2 ,...,n i ,...,n M (t)\ 2 (13) 

ni =0,Vi^/3,7 

where the sum over rii does not include np or n 7 , as they 
are held fixed and particle conservation requires a fixed 
total number of particles. P(n^,n 7 ,t) shows the proba- 
bility of finding np particles in (3th well simultaneously 
with finding n 7 particles in the 7th well. 

A. Large entangled number states in four wells 




FIG. 10: (color online). Ground and entangled number state 
for 128 particles in four wells: (a) Ground state and (b) en- 
tangled number state at t = 2 AO ms for 128 particles in four 
wells, U OMEr, J 0.04e- a E R with a 0.175. The 
time shown in (b) represents the earliest time when maximum 
probability occurs for an entangled number state. Shown is 
the joint probability function P(ni, 712) versus the number of 
particles in wells 1 and 2 on the x and y axes, respectively. 

The results of our investigation of the creation of en- 
tangled number states in four wells with a large number 
of particles in each well are provided below. The regime 



FIG. 11: (color online). Varying the tunneling parameter, J, 
for four wells: (a) Ground state and (b) entangled number 
state for 128 particles in four wells with U = 0. 01 Er and 
J 0Me~ a E R with a 0.175 and t 2.40 ms, a 1.175 
and t = 3.65 ms, and a = 2.175 and t — 5.70 ms in rows 1, 
2, and 3. This shows that initial squeezing can be used to 
control the final entanglement as have been shown for three 
wells in Fig. 7 and double wells in Ref. [24] 



we explore in four wells is similar to the Kasevich [23] ex- 
periments and we report results using experimental pa- 
rameters relevant to their work for a condensate made of 
87 Rb with A = 840 nm and a sc = 5.8 nm. We assume an 
intial phase offset of tt between the wells. Fig. 10 shows 
the joint probability function for (a) the ground state at 
t = and (b) entangled number state at time t = 2.40 ms 
for 128 particles in four wells. In this figure and all follow- 
ing the time depicted represents the earliest time when 
maximum probability for an entangled state occurs. The 
ground state has a Gaussian distribution centered around 
the state where the particles are equally distributed be- 
tween the wells. The joint probability function for the 
entangled number state has two peaks one where there 
are 57 particles in the (3 = 1 well and 7 particles in the 
7 = 2 well and another peak of equal probability where 
there are 7 particles in the (3 = 1 well and 57 particles 
in the 7 = 2 well. Besides this highest probability state, 
there are other non- vanishing coefficients that are smaller 
and distributed around this peak. 

Fig. 11 shows the effect of varying J, the tunneling pa- 
rameter. Column (a) shows the joint probability function 
for the ground state at t = and column (b) shows an 
entangled number state at t = 2.40, 3.65, and 5.70 ms in 
rows 1, 2, and 3, respectively. The tunneling parameter 
is set to J = 0.04e- a E R with a = 0.175, a = 1.175, and 
a = 2.175 for rows 1,2, and 3 respectively. Increasing a 
corresponds to decreasing tunnelling and correlates with 
increased time to evolve into the entangled number state. 
Increased a also means starting with a state with larger 
number squeezing, and thus we see in panels (a), (b) and 
(c) that the extremity of the final cat state is determined 
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by initial squeezing. This same behavior was shown with 
double well [24] , and with three wells in the previous sec- 
tion, and thus we show that it is a general characteristic 
of multi-well cat states. We should emphasize the find- 
ing that to obtain a more extreme cat state, the initial 
barriers have to be low. 

Fig. 12 illustrates the effect of varying the number of 
particles. Column (a) shows the ground state at t = 
and column (b) shows the entangled number state for 256 
particles at t = 1.74 ms, 384 particles at t = 1.47 ms, and 
512 particles at t = 1.30 ms. Increasing the particle num- 
ber raises the strength of the interaction term in the BH 
Hamiltonian and results in decreased tunneling, explain- 
ing the increased squeezing of the ground and entangled 
number states with increased particle number. 



(a) 



(a) 

£ 

1 
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(c) 



. 25Z ^ 6 19£ 



(d) 1 



2 <f) 



i i 



on m I' 



FIG. 12: (color online). Varying number of particles for four 
wells: (a) Ground state and (b) entangled number state with 
U = OMEr and J = 0.04e~ a E R with a = 0.175 for 256 
particles at t — 1.74 ms, 384 particles at t = 1.47 ms, and 512 
particles at t = 1.30 ms in four wells for rows 1,2, and 3. For 
the same barrier height, higher number of particles give rise 
to less extreme superpositions. 

Our simulations indicate that entangled number states 
evolve through natural time evolution of the BEC with 
an initial phase offset of ir between the wells even with a 
large number of particles in each well. The extremity of 
the entangled number state can be controlled by varying 
the tunneling parameter or the number of particles. 



B. Entangled number states in eight wells 

We present the dynamic creation of entangled number 
states in eight wells with just a few particles per well; this 
regime is similar to the regime studied by Greiner [22] 
and we report our results using experimental parameters 
relevant to their work; 87 Rb, A = 852 nm, a sc = 5.8 
nm. Fig. 13 shows (a) ground state and (b) the time 
evolved entangled number state for 24 particles in eight 
wells with J 0Me~ a E R with a 0.175 and U 



(b) 



FIG. 13: (color online). 24 particles in eight wells: (a) Ground 
state and (b) entangled number state at t = 2.79 ms for 24 
particles in eight wells with U = 0.0 IEr and J = 0.04e~ a Er 
with a = 0.175. 



O.Oli^R with an initial phase offset of tt between the wells. 
Our investigations show that entangled number states 
can also be realized in a larger number of wells with a 
small number of particles in each well. 



C. Analysis of entangled number states: a new 
type of entangled state 

The extreme entangled state for a multiwell BEC takes 
the form of Eq. 2. And the less extreme entangled states 
in multiple wells can have the approximate form 



\^} {N,M} non - extr 



(\N - jn,n, ...,n) 



+ \n,N -jn,n,...,n) 

+ ... + |n,...,7V-jn))(14) 

where n << N and j = M — 1. The four 
well cat state shown in Fig. 8(b) is an ex- 
treme cat state in the approximate form = 
^ (|16, 0, 0, 0} + |0, 16, 0, 0} + |0, 0, 16, 0) + |0, 0, 0, 16)). 



An 

state 



example 
would 
1,1) + 



of a 
be of 
|1,13,1,1)- 



cat 



less extreme 
the form |^) = 
-|1,1,13,1) + |1,1,1,13)). 



Thus, one would expect a similar form for the large 
entangled states described in this section, however this 
is not what we found, as is shown below. 

The states in Fig. 10 with the highest probability at 
t = and t = 2.40 ms are shown in the first two entries 
of Table I and indicate there are 32 particles in each of 
the wells in the ground state and 57 particles simultane- 
ously in all four wells in the entangled number state. The 
highest probability states in Fig. 11 with varying J are 
shown in Table I. Table I shows the ground states with 
the particles equally distributed between the wells and 
the entangled states with 57, 48, and 42 particles in all 
the wells. The highest probability states for Fig. 12 with 
varying TV are shown in Table II, and Table III shows 
the highest probability states for the ground and entan- 
gled number states for the eight well system shown in 
Fig. 13. We would like to emphasize that the listed state 
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is the highest coefficient state, and there are other lower 
coefficient states around this. 




FIG. 14: (color online). Evolution of entangled state for 128 
particles in four wells, shown in the joint probability distribu- 
tion space of Eq. 13 for sites 1 and 2. Shown are (a) at t=0, 
the ground state, (b) at 1—1.22 ms, the ground state broad- 
ens, (c) at t=1.71 ms, (d) at t=1.84 ms the state is beginning 
to split just as it does in a double well and triple well, (e) at 
t=1.96 ms, and (f) at 1—2A2 ms, the final entangled state is 
formed. The evolution here further illustrates that the under- 
lying physical mechanism for the creation of cat states is the 
same in a double well and multiple wells. 

Based on the empirical evidence of the entangled num- 
ber states shown in this section we find the entangled 
number states for large N/M and large M tend towards 
the approximate form 




where n < N and there are ^ — n particles simul- 
taneously in all M wells. We say 'approximate form' 
because although the highest probability state has this 
form, there are many non-vanishing coefficients as well 
with smaller amplitudes, and spread out around the high- 
est probability peak. The initial parameters of the sys- 
tem and dynamics determines the relative amplitudes 



TABLE I: Highest probability states for 128 particles in four 
wells with U — O.OIEr and varying J, where J — 0.04e _Q! Er. 
These states correspond to Fig. 11. It should be emphasized 
that all the Tables in this article show only the highest state; 
however, there are other non- vanishing coefficients as well 
as evident from Fig. 11. Those coefficients are smaller and 
spread out around the highest probability peak. 

t (ms) a highest probability state 
0.0 0.175 |32,32,32,32) 
2.4 0.175 ^(|57, 7, 57, 7) + |7, 57, 7, 57)) 
0.0 1.175 |32,32,32,32) 
3.65 1.175 ^=(|48,16,48, 16) + |16,48,16, 48)) 
0.0 2.175 |32,32,32,32) 

5.7 2.175 4-(|42,22,42,22) + |22, 42, 22, 42)) 



TABLE II: Highest probability states in four wells with J = 
0.04e-° 175 E R , U = 0ME R and varying N. These states 
correspond to Fig. 12. 

t (ms) N highest probability state 

0.0 256 |64,64,64,64) 

1.74 256 t=(|102, 26, 102, 26) + |26, 102, 26, 102)) 

0.0 384 |96,96,96,96) 

1.47 384 ^(|143,49, 143,49) + |143,49,143,49)) 

0.0 512 |128,128,128,128) 

1.30 512 -4-Q183, 73, 183, 73) + |73, 183, 73, 183)) 



of the non-vanishing coefficients for the final entangled 
states generated in our model. 

To our knowledge at the time of writing, number entan- 
gled states in the form of Eq. 15 have not been discussed 
elsewhere and represent a new type of entanglement. 
Well known entangled states include Greene, Home, Zeil- 
nger states [15] shown in Eq. 1 which for two wells and 
TV particles is an extreme entangled state of the form, 
(l/v / 2)(|7V,0) + |0,7V)), also known as NOON states. An- 
other type of state is the W-state [46] , which for TV qubits 
has the form, (1/ V / 7V)(|10..00) + |01..00) + .. + |00..01)). 
Because of the 'twin-like' nature, we refer to the num- 
ber entangled states of Eq. 15 as Siamese states. For the 
case of qubits, W state was found to be more robust than 
GHZ state. For multi-dimensional quDits, the Siamese 
states introduced here may have different properties than 
other forms of number entangled states. We have only 
shown here a way of generating such novel states with 
cold atoms in a lattice without yet knowing their prop- 
erties. 

Tables IV and V illustrate the relationship between 
particle number and interaction energy and these new 
types of entangled states. Table IV shows that when 
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TABLE III: Highest probability states in eight wells with 
J = 0.04e~°- 175 E R , U = O.OIEr. These states correspond 
to Fig. 13. 

t (ms) highest probability state 
0.0 |3,3,3,3,3,3,3,3) 
2.79 -U|6, 0, 6, 0, 6, 0, 6, 0) + |0, 6, 0, 6, 0, 6, 0, 6)) 



we increase the number of particles, the Fock state num- 
ber distribution in the final evolved state slowly turn into 
twin-like states. Same observation can also be made from 
the increase of interaction parameter U. Since particle 
number increment is effectively an increase in effective in- 
teraction, these results complement each other. Since in 
the Bose-Hubbard model the properties are determined 
by the ratio of U/t, we can say that starting with a 
small value of tunneling or a high barrier (highly squeezed 
state) is a prerequisite for obtaining these states. This ac- 
tually makes the distinction very clear for the four well 
results, since we have shown for double well and three 
well, that to get toward an extreme cat state limit we 
need to start in the opposite limit of low barrier (gaus- 
sian state). 

Fig. 14 shows evolution of a four well state, for a twin- 
like state. In the graphical representation of two-site 
joint probability function, a twin-like state and a four- 
legged cat state cannot be differentiated. For that we 
look at the coefficients of different Fock states to verify 
that this is a twin-like state, with the highest probability 
state being (|57, 7, 57, 7) + |7, 57, 7, 57)), with other non- 
vanishing coefficients that are smaller. Fig. 14 presents 
the general picture of how a four well state evolves - 
(a) shows the phase imprinted ground state, (b) and (c) 
shows the broadening of the Gaussian distribution, (d) 
and (e) shows the process of splitting into superposi- 
tions, and (f) shows the final evolved state. This en- 
tangled number state example is for 128 particles, and 
U = 0.01E R , J = 0Me- a E R with a = 0.175. Figure 14 
can be compared with Figs. 4 and 6 for double well and 
triple well respectively. We see that the dynamics of en- 
tangled state generation in a double well and multiple 
well follow similar properties - a gaussian ground state 
once phase imprinted to an unstable equilibrium, evolves 
first by broadening the gaussian, and then slowly split- 
ting symmetrically to a macroscopic superposition state. 

The detailed mechanism whereby the entangled states 
of Eq. 2 cross over to the new type of entangled states, 
with alternating populations as observed here, is not fully 
understood. Nor it it obvious why other types of entan- 
glement pairings and types of alternation do not seem 
to occur. However, the most likely explanation of the 
cross-over effect itself will simply involve the fact that 
the entangled states of Eq. 2 have the maximal possible 
energies for the system: the proposed types of phase en- 
gineering explored here simply do not lead, as TV and U 



increase, to energies high enough to allow simple time 
evolution into entangled states of such high energy. In 
terms of the underlying classical phase space dynamics, 
we can say that the characteristics of the unstable equi- 
librium (0,0,0,7r,7r,7r) (in the conjugate variables ni2, n23, 
^34, 012, 023, 034) changes as U and N is increased. For 
the double well case we know that U N/ J > 1 makes the 
7r-phase an unstable equilibrium, and stable otherwise. 
Similarly, here the combination of values of U and N 
changes the scenario and beyond a certain value of UN, 
a complicated phase space dynamics constrains the mo- 
tion in phase space in such a way that only certain region 
is accessible. Since we have presented quantum dynam- 
ics of multi-well in this paper based on intuition from its 
semiclassical aspects without a detailed study of dynam- 
ics in the underlying classical phase space as was done 
for the double well [24], we cannot quantify our findings. 

TABLE IV: Highest probability states in four wells with J = 
0.04e~ a E R with a = 0.0, U = 0.01E R and varying N. 

t(ms) N highest probability state 

7.69 20 ^(|20, 0,0,0) + |0,20,0,0) + 

|0,0,20,0) + |0,0,0,20)) 
6.39 24 ^(|23,0, 1,0) + |0,23,0, 1) + 

|1,0,23,0) + |0,1,0,23)) 
6.24 28 ^(|26,0,2,0) + |0,26,0,2) + 

|2,0,26,0) + |0,2,0,26)) 
5.38 32 ^=(|27, 0,5,0) + |0,27,0,5) + 

|5,0,27, 0) + |0,5,0,27)) 
2.91 64 ^(|32, 0, 32, 0) + |0, 32, 0, 32)) 
2.15 128 ^=(|59,5,59,5) + |5, 59, 5, 59)) 



VI. CONCLUSIONS 

We have demonstrated phase engineering schemes for 
the generation of entangled number states in multiple 
wells. These states represent a superposition of particles 
in a large number of spatial locations - a multi-positional 
generalization of double well NOON states. It is shown 
that entangled number states can be evolved from phase 
imprinted ground states of BEC in multiple wells, and the 
number of particles participating in the entanglement can 
be controlled by varying the height of the barrier between 
the wells, controlling the rate of ramping, and/or the 
number of particles. The less extreme entangled number 
states, which we show how to create in a controlled way, 
represent macroscopic superposition states that are more 
robust to experimental conditions and particle loss. We 
demonstrated our scheme with small number of particles 
in a small lattice, large number of particles per well, as 
well as with a large number of wells. Our investigations 
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TABLE V: Highest probability states in four wells with J — 
0Me~ a E R with a = 0.0, N = 32 and varying U with U = 
u(0MEr). 

t(ms) u highest probability state 

5.70 0.15 ^(|32, 0,0,0) + |0,32,0,0) + 
|0,0,32,0) + |0,0,0,32)) 
5.38 0.25 7|(|27, 0, 5, 0) + |0, 27, 0, 5) + 
|5,0,27, 0) + |0,5,0,27)) 
2.91 0.50 ^(|18, 0,14,0) + |0,18,0,14) + 
|14, 0,18,0) + |0,14,0,18)) 
2.47 0.75 ^(|17,1,13,1) + |1,17,1,13) + 
|13, 1,17, 1) + |1,13,1,17)) 
2.08 1.00 ^=(|16,2,12,2) + |2, 16,2, 12) + 
|12,2,16,2) + |2, 12, 2, 16)) 
1.66 1.25 tj(|14, 2, 14, 2) + |2, 14, 2, 14)) 



equilibrium splits the wavepacket symmetrically to cre- 
ate these states. We briefly describe our earlier work on 
double well [24, 25] to demonstrate this in the semiclassi- 
cal phase space, and to motivate and generalize the study 
to larger number of wells and lattices. The significance 
of the phase imprint values of 27r/3 for three wells, and 7r 
for four wells is explained this way. For the initial phase 
difference between the neighboring wells required to cre- 
ate these states, we presented a novel series of formulae 
that is valid for any number of wells, even or odd. 

The creation, characterization, and applications of 
multidimensional/multipositional Schrodinger cat states 
of atoms discussed in this article remain largely unex- 
plored experimentally, and the theoretical ramifications 
of such states, should they be easily produced, are just 
emerging. 



Acknowledgments 



for four wells revealed a new type of number entangled 
state where a twin like state is formed, which we refer to 
as Siamese states. We demonstrated a relationship be- 
tween the particle number and interaction strength and 
the formation of these new types of entangled states. 

The physical mechanism by which these states are gen- 
erated can be understood in terms of the underlying semi- 
classical phase space. Ground state put on an unstable 



This work was supported by NSF grant PHY-0140091 
and PHY-07-03278 and MAL gratefully acknowledges 
support from the DOE computational science gradu- 
ate fellowship program grant DE-FG02-97ER25308 and 
used resources of the National Energy Research Scientific 
Computing Center, which is supported by the Office of 
Science of the U.S. Department of Energy under Contract 
No. DE-AC03-76SF00098. 



[1] A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 
777 (1935). 

[2] J. S. Bell, Speakable and Unspeakable in Quantum Me- 
chanics (Cambridge Univ. Press, Cambridge, 1987). 

[3] R. Horodecki, P. Horodecki, M. Horodecki, and K. 
Horodecki,, Rev. Mod. Phys. 81, 865 (2009). 

[4] S. F. Huelga et a/., Phys. Rev. Lett. 79, 3865 (1997); V. 
Giovannetti, S. Lloyd, L. Maccone, Phys. Rev. Lett. 96 
010401 (2006); 

[5] L. Pezze and A. Smerzi, Phys. Rev. Lett. 102 100401 
(2009). 

[6] C. Monroe, Nature 416, 238 (2002). 
[7] Z. Zhao et a/., Nature 430, 54 (2004). 
[8] H. Haffner et a/., Nature 438, 639 (2005). 
[9] J. A. Jones et al, Science 324, 1166 (2009). 
[10] O. Mandel, M. Greiner, A. Widera, T. Rom, T. W. Han- 
sen, I. Bloch, Nature 425, 937 (2003). 
[11] J. Esteve, C. Gross, A. Weller, S. Giovanazzi, and M. K. 

Oberthaler, Nature 455, 1216 (2008). 
[12] W. Marshall, C. Simon, R. Penrose, and D. 

Bouwmeester, Phys. Rev. Lett. 91, 130401 (2003). 
[13] O Romero-Isart, M. L. Juan, R Quidant, and J. I. Cirac, 

New J. Phys. 12, 033015 (2010). 
[14] J. S. Bell, Physics 1, 195 (1965). 

[15] D. M. Greenberger, M. A. Home, A. Shimony, and A. 
Zeilinger, Am. J. Phys. 58, 1131 (1990). 



[16] E. Schrodinger, Naturwissenschaften 23, 807 (1935). 
[17] J. R. Friedman et a/., Nature 406, 43 (2000). 
[18] C. H. van der Wal et al, Science 290, 773 (2000). 
[19] T. Morimae, Phys. Rev. A 81, 010101(R) (2010). 
[20] N. J. Cerf, S. Massar, and S. Pironio, Phys. Rev. Lett. 

89, 080402 (2002). 
[21] O. Morsch, M. Oberthaler, Rev. Mod. Phys. 78, 179 

(2006); M. Leweinstein, A. Sanpera, V. Ahufinger, B. 

Damski, A. Sen, and U. Sen, Advances in Physics 56, 

243 (2007). 

[22] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and 

I. Bloch, Nature 415, 39 (2002). 
[23] C. Orzel, A. K. Tuchman, M. L. Fensclau, M. Yasuda, 

and M. A. Kasevich, Science 291, 2386 (2001). 
[24] K. W. Mahmud, H. Perry, and W. P. Reinhardt, Phys. 

Rev. A 71,023615 (2005); K. W. Mahmud, H. Perry, and 

W. P. Reinhardt, J. Phys. B 36, L265 (2003). 
[25] K. W. Mahmud, Ph.D. Thesis, University of Washington 

(2004). 

[26] L. Amico, A. Osterloh, and F. Cataliotti, Phys. Rev. 
Lett 95, 063201 (2005); K. Henderson, C. Ryu, C. Mac- 
Cormick and M. G. Boshier, New J. Phys. 11, 043030 
(2009). 

[27] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. 

S. Fisher, Phys. Rev. B 40, 546 (1989). 
[28] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 



13 



P. Zoller, Phys. Rev. Lett 81, 3108 (1998). 
[29] M. A. E. Leung and W. P. Reinhardt, Comp Phys Comm, 

177 (4), 348 (2007). 
[30] J. Dunningham and D. Hallwood, Phys. Rev. A 74, 

023601 (2006); J. A. Dunningham, K. Burnett, R. Roth, 

and W. D. Phillips, New J. Phys. 8, 182 (2006). 
[31] A. Nunnenkamp, A. M. Rey, and K. Burnett, Phys. Rev. 

A 77, 023622 (2008). 
[32] P. W. Anderson, Lectures on The Many- Body Problem, 

Vol 2, Pages 113-135, Academic Press, 1964. 
[33] A. Smerzi, S. Fantoni, S. Giovanazzi, and S. R. Shenoy, 

Phys. Rev. Lett. 79, 4950 (1997); S. Raghavan et al, 

Phys. Rev. A 59, 620 (1999). 
[34] M. Albiez, R. Gati, J. Foiling, S. Hunsmann, M. Cris- 

tiani, M. K. Oberthaler, Phys. Rev. Lett. 95, 010402 

(2005); 

[35] K. Husimi, Proc. Physico-Math. Soc. Japan 22, 264 
(1940); H. Lee, Phys. Rep. 259, 147 (1995). 

[36] J. Denschlag et al, Science 287, 97 (2000). 

[37] A. Micheli, D. Jaksch, J. I. Cirac, and P. Zoller, Phys. 
Rev. A 67, 013607 (2003). 

[38] F. Piazza, L. Pezze, and A. Smerzi, Phys. Rev. A 78, 



051601 (2008). 

[39] A variety of array configurations have been studied, in- 
cluding linear and circular arrays. The trimeric linear 
array with differing a has been studied by: P. Buon- 
sante, R. Franzosi, V. Penna, Phys. Rev. Lett 90, 050404 
(2003). 

[40] R. Franzosi and V. Penna, Phys. Rev. E 67, 046227 
(2003). 

[41] S. Mossmann and C. Jung, Phys. Rev. A 74, 033601 
(2006). 

[42] F. Trimborn, D. Witthaut, and H. J. Korsch, Phys. Rev. 

A 79, 013608 (2009). 
[43] T. F. Viscondi, K. Furuya, M. C. de Oliveira, EPL 90, 

10014 (2010). 

[44] A. Polkovnikov, Phys. Rev. A 68, 033609 (2003); A. 
Polkovnikov, S. Sachdev, S. M. Girvin, Phys. Rev. A 66, 
053607 (2002). 

[45] K. W. Mahmud, M. A. Leung, and W. P. Reinhardt, 

arXiv : cond-mat/0403002. 
[46] W. Dur, G. Vidal and J. I. Cirac, Phys. Rev. A 62, 

062314 (2000). 



